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In estimating annual construction budgets and work forces required 
for a large entity, such as an operating company, it is desirable to know 
the total annual number of cable pairs to be added at the main frames 
in all of the feeder routes in the entity. These pairs are the result of 
many independent relief decisions made in the individual feeder routes. 
A model is presented which relates the total annual number of cable 
pairs added at the main frames to aggregate relief timing and sizing 
design parameters, an aggregate demand forecast and the total as- 
signed and available pairs at the main frames in all of the routes in the 
entity. Because the model does not require individual route data or 
specific relief projects, it can be used as a "top-down" check on the 
"bottom-up" requirements determined by aggregating a list of iden- 
tified relief projects. It is also useful for estimating cable pair re- 
quirements when specific relief projects are not yet identified. 

I. INTRODUCTION 

Growth in the number of telephone subscribers requires the periodic 
addition of cable pairs in the feeder network. 1 The feeder network 
comprises the large backbone cables which funnel cable pairs from the 
distribution network back toward the local wire center. The feeder cables 
are terminated in the wire center on a main distributing frame which 
serves as the interface between the loop network and the switching 
equipment. 

At the main frame, the feeder cables are grouped into separate feeder 
routes which serve disjoint geographical areas within the boundaries of 
the wire center. Each route at the main frame is composed of a number 
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of available pairs, some of which are already assigned to existing sub- 
scribers served by the feeder route. The ratio of assigned pairs to avail- 
able pairs at the main frame is called the main frame fill of the feeder 
route. The fill at other points along the route is defined in a similar 
manner. 

When the fills in a portion of the feeder route become too high, new 
cable must be added. Most relief cables extend cable pairs already ter- 
minated on the main frame further out into the route. If the main frame 
fill of the route is also high, the new cable will be terminated on the main 
frame. This relief process is a complex one, involving both economic 2 and 
physical 3 considerations in the particular feeder route. 

Individual relief timing 4 and sizing 2 decisions are made throughout 
the year by engineers based in local district engineering offices. A typical 
district has responsibility for about 100 main feeder routes. Districts are 
combined into administrative entities called areas which contain about 
4 districts and 400 feeder routes. A typical operating company contains 
about 4 areas and 1600 feeder routes. 

The aim of this paper is to present a model which has been derived 
to estimate the total number of cable pairs added on the main frames 
during a year in a large entity such as an area or company. The number 
of these pairs has historically been related to the total amount of feeder 
cable used in the entire feeder network during the year, and hence are 
an important component of an entity's annual construction budget and 
work force estimates. The model uses aggregate timing and sizing design 
parameters, an aggregate demand forecast and known aggregate main 
frame data such as total available pairs and total assigned pairs in the 
entity. 

The model can be used to estimate future cable pair requirements for 
a large entity without identifying individual feeder routes or individual 
relief jobs. Hence the model is useful for estimating future cable pair 
requirements when specific relief projects are not yet identified, and for 
providing a simple "top-down" check on the "bottom-up" requirements 
obtained by aggregating a "market list" of identified relief projects. This 
check is useful because of the large number of feeder relief projects re- 
quired each year in a large entity such as an area or company. 

II. SUMMARY 

A model is derived to estimate the total increase in available pairs at 
the main frames next year in all of the routes in a large entity. The model 
can be applied recursively to estimate the total increase in available pairs 
for any future year. 

The problem of concern is formulated in Section HI. Known aggregate 
main frame data are defined for an entity in terms of analogous data from 
the individual feeder routes. The basic approach used to derive the model 
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is outlined in Section 3.3 and involves modeling route relief at the main 
frame as a binomially distributed random variable where success con- 
notes terminating a new cable on the main frame next year. Most of the 
modeling effort involves calculating the probability of relief (success) 
next year. Assumptions used to derive the model are conveniently 
summarized by defining an "analog entity" in Section 3.4. 

The model is derived for an analog entity in three stages in Section 
IV. The initial result assumes that the fill at relief — the main frame fill 
in the route when a new cable is terminated — and cable size terminated 
at next relief will be the same as the fill at relief and cable size terminated 
at last relief in all of the routes in the analog entity. Adjustments to this 
steady state result are then derived in Sections 4.2 and 4.3 which enable 
the model to account for both changes in fill at relief and changes in cable 
size terminated. 

The general model for estimating the increase in available pairs (AP) 
next year in an entity is 

SP - j + (1 - B)P (l - j\ + AT AS 

where 

A- W 



P--NS 
2 



is an estimate of the average fill at last relief for the entity. The estimates 
are expressed in terms of known aggregate main frame parameters for 
the entity which are either observed at the beginning of next year or 
forecast for next year: 

N = number of routes 

W = total assigned pairs 

P = total available pairs 

S = average cable size terminated 

G = growth in assigned pairs next year 

F = average fill at next relief 

AS = change in average cable size terminated next year 

The parameter depends on the impedance of the entity to a change in 
fill at relief. 

Confidence intervals are calculated for the model in Section 4.5 as a 
function of the number of routes in the entity for the idealized case of 
identical routes. The resulting error bounds represent the best possible 
model performance in actual entities. 

Finally, a validation study is described in Section V where the per- 
formance of the model is evaluated in actual entities. This study also 
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provided the means for empirically selecting a universal value of 0.6 for 
the impedance parameter 0. An error bound of 12 percent at the 90 
percent confidence level was determined for the = 0.6 model when 
applied to a typical company-sized entity of 1600 routes. The corre- 
sponding error bound at the 50 percent confidence level, or probable 
error, was 5 percent. 

III. MODEL FORMULATION 

A model for estimating the total increase in available pairs at the main 
frames in a given entity next year will be derived. This parameter is an 
indicator of the total feeder cable requirements in the entity next year. 
The goal of the modeling effort is to estimate this parameter using known 
aggregate main frame data and to determine the statistical accuracy of 
the estimate when used in actual entities. 

Known aggregate main frame data for an entity include: 

• number of feeder routes 

• total assigned pairs 

• total available pairs 

• average cable size terminated 

• growth forecast 

• average fill at next relief 

• change in average cable size terminated 

These data are either aggregated annually from the feeder routes which 
compose the entity, or are forecast directly for the entity. The above 
parameters will be formally defined in terms of analogous route pa- 
rameters in Section 3.2. 

3. 1 Basic route parameters 

A feeder route is described at the main frame by the following set of 
basic route parameters (denoted by lower-case letters) which are either 
observed or forecast at the beginning of next year (see Fig. 1): . 

w = assigned pairs \ 

p = available pairs I observed at the 

s = cable size terminated j beginning of next year 

a = fill at last relief ) 

g = growth in assigned pairs next year 

/ = fill at next relief 

The ratio w/p is the main frame fill at the beginning of next year. The 
last cable terminated at the main frame contained s pairs and was placed 
in service when the main frame fill, or fill at last relief, was a. A new cable 
will be terminated in the route when the main frame fill reaches the fill 
at next relief, /. Initially, we will assume that the new cable contains the 
same number of pairs as the last cable terminated (s). This assumption 
will be relaxed in Section 4.3. 
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Fig. 1 — Feeder route parameters. 



3.2 Known aggregate main frame data 

Known aggregate main frame data (denoted by upper-case letters) 
compiled from the individual feeder routes at the beginning of next year 
are 

N = number of feeder routes 
W = Y. w j ~ total assigned pairs 
P = YlPj = total available pairs 

S = — Y.Sj = average cable size terminated 

The parameters W and P are simply aggregate facility data which have 
historically been collected to monitor the usage of the feeder network 
at the main frames. The ratio W/P is the main frame fill of the entity at 
the beginning of next year. The average cable size terminated (S) is an 
aggregate sizing design parameter. 

The following parameters are forecast for the entity for next year: 

G = J2ij = growth in assigned pairs next year 

F = — Y,Pjfj = average fill at next relief. 

G is a standard demand forecast and F is an aggregate timing design 
parameter. These parameters are typically forecast directly for the entity 
rather than aggregated from the forecasts of the analogous parameters 
in each route. 

The definition of N, W, P and G are intuitively aggregates of the 
analogous feeder route parameters. The definition of the aggregate de- 
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sign parameters S and F will become clear during the derivation of the 
model (Section IV). Note that there is no known aggregate parameter 
which is analogous to the fill at last relief (a) in a single route (see Section 
3.1). This unknown aggregate parameter will be introduced with the 
definition of the analog entity in Section 3.4. 

The basic approach used to derive a model for estimating the total 
increase in available pairs next year in an entity is described in the next 
section. 

3.3 Basic approach 

The increase in available pairs (Ap) next year in a route is either s or 
zero, depending on whether a new cable is terminated next year or not. 
Hence Ap can be expressed as x s where x is the relief indicator for the 
route next year, defined as 



X = 



1 if a new cable is 

terminated next year 
otherwise 



The increase in available pairs next year in an entity can be expressed 
as 

AP=£Ap; 

= ZxjSj (1) 

The minimum value of AP = corresponds to the unlikely case of no new 
cable being terminated next year (%>• =0:j = l,iV), while the maximum 
value of AP = NS corresponds to the equally unlikely case of a new cable 
being terminated in every route next year (xj = 1 : j = 1,N). 

The relief indicator (x) for a route next year can be modeled as a bi- 
nomially distributed random variable where success corresponds to a 
new cable being terminated next year (x — 1). Treating the cable size (s,-) 
as a known parameter, the expected increase in available pairs next year 
in the entity can be expressed as 

AP=£X,sy (2) 

where X y = xj is the probability of relief next year for route /. 

Intuitively, the probability of relief next year should depend on how 
far the route is through its relief cycle at the beginning of next year. The 
fraction through the relief cycle (r) decomposes the relief cycle (t) into 
two components; the time since last relief (rt) and the time to next relief 
[(1 — r)t]. The route will be relieved next year if the time to next relief 
(observed at the beginning of next year) is one year or less. If r is modeled 
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as a random variable, then the probability of relief next year, 

= P (^l-i) (3) 

depends on the tail of the distribution of the fraction through the relief 
cycle at the beginning of next year. 

The route is said to be in the steady state if the fill at next relief is the 
same as the fill at last relief (/ = a). In this case (see Fig. 2), the relief cycle 
of the route — the number of years between last and next relief — is 



sa 

t ss - 



(4) 



and the fraction through the relief cycle at the beginning of next year 
is 



pSS — 



w — (p — s)a 



sa 



(5) 



If the fill at next relief (/) forecast at the beginning of next year is not 
equal to the fill at last relief (a), then the relief cycle becomes 



t = 



p(f -a) + sa 

g 



(6) 



If / > a (see Fig. 3), then the relief cycle is increased and the fraction 
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through the relief cycle is decreased such that their product, rt = time 
since last relief, remains constant. Because the time since last relief is 
also equal to r ss t ss , we have rt = r ss t ss or 



r = r» 



(7) 



If r sa is modeled as a random variable, then the probability of relief 
next year is 



*-f.(r*l-i) 

= P (r ss > ^ 

V t ss t ss / 



(8) 



Hence r ss — the fraction through the relief cycle if the route is in the 
steady state (/ = a) — is the basic random variable which will be used to 
derive the model. 

The distribution of r ss determines the distribution of w — the assigned 
pairs in the route at the beginning of next year — which can be expressed 
in terms of r ss by solving eq. (5) for w, 



w = ap + as(l — r ss ) 



(9) 



This is the only basic route parameter defined in Section 3.1 which is 
modeled as a random variable. All other parameters (p,s,a,g,f) are as- 




YEARS 

Fig. 3 — Feeder route with a change in fill at relief (/ > a). 
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sumed to be deterministic. Hence the length of the relief cycle [t ss and 
t, see eqs. (4) and (6)] are deterministic, but the fraction through the 
relief cycle [r ss and r, see eqs. (5) and (7)] are random variables. 

3.4 Model assumptions — the analog entity 

Several assumptions are required to derive a model for the total in- 
crease in available pairs (AP) next year in an entity in terms of the known 
aggregate main frame data (N,W,P,S,G,F) defined in Section 3.2. These 
assumptions are needed because the distribution of the aggregate data 
over the routes in the entity (wj, pj, Sj, a,, gj, fj : j = 1,N) is in general not 
known. The assumptions are: 

Al: All of the routes in the entity are assumed to have the same fill at 
last relief (aj = A:j = 1,N) where A is an unknown parameter that 
must be estimated in terms of known aggregate data. 

A2: If the routes in the entity are in the steady state condition (fj = A 
: ; = 1,N), then the fractions through the relief cycle of the routes 
(r] s : j = 1,N) are independent and identically distributed uniform 
random variables on the unit interval. In other words, if a route is 
selected at random from routes in the steady state, then we are 
equally likely to observe any value of r ss between (just after relief) 
and 1 (just before relief). 

An entity satisfying these two conditions will be referred to as an an- 
alog entity. The model will be derived for an analog entity in the next 
section. The performance of this model in actual entities will be studied 
in Section V. 

IV. AVAILABLE PAIR INCREASE MODEL 

A model for estimating the increase in available pairs next year in an 
analog entity is derived. This model can be applied recursively to obtain 
an estimate of the increase in available pairs for any future year. 

4.1 Steady state model 

The expected increase in available pairs in the analog entity next year 
[eq. (2)] is 

AP=EX,s, 

where Sj is the cable size terminated and Xj is the probability of relief 
next year in route j. If the routes in the analog entity are in the steady 
state condition (i.e., fj = A:j = 1,N), then tj = tf [see eqs. (4) and (6)] 
and the probability of relief next year [eq. (8)] is 



*-*(**%-$ 

-*(**!-£) 
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Since r s / is assumed to be uniformly distributed on the unit interval 
(assumption A2), we have 



*»-■ 



1 if tf<l 

1 (10) 

— if tr>i 

+ SS J 



The physical significance of tf < 1 is that the last cable terminated (sj) 
was sized so small that the route requires another relief cable in less than 
one year. This is not a normal engineering practice and it will not be 
considered in the model. 

The expected available pair increase in route j for the case tf > 1 
is 



Apj = XjSj 

_ 1L 

4.SS 



(11) 



.a 

A 

This result follows from eqs. (10) and (4) and assumption Al. Hence if 
Sj pairs are terminated every t f years, then Sj/tf pairs will be terminated 
on average next year. Equivalently, if cables were sized for only one year's 
growth, then gj/A pairs would be terminated next year. 

The expected increase in available pairs in the analog entity next year 
is obtained by summing eq. (11) over all routes, 

A 

= f (12) 

A 

where G is the entity growth forecast for next year. In the following 
sections, we derive adjustments to this simple steady state model by 
relaxing some of the restrictive assumptions leading to this result. 

4.2 Change In fill at relief 

The steady state assumption will now be relaxed by allowing the 
forecast fill at next relief in route ; (fj) to be different from the fill at last 
relief (A). Such a change could take place if locally selected fills at relief 
were used in place of a constant fill at relief policy represented by A. 
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4.2. 1 Instantaneous change in fill at relief 

If the fill at next relief is changed from A to / ; at the beginning of 
next year, then the probability of relief next year in route j [eq. (8)] is 

1 \ ' ty ty) 

(0 if tj > ty + 1 

VI if tj < 1 

This result follows because r] s is uniformly distributed on the unit in- 
terval by assumption A2. If fj » A such that a year or more is added to 
the relief cycle (case tj > ty + 1), then the route will not be relieved next 
year (X; = 0). On the other hand, if fj « A such that the relief cycle is 
reduced to less than one year (case tj < 1), then the route will be relieved 
next year (X ; = 1). 

The expected increase in available pairs in route j for the case 1 < tj 

< ty + i is 

Apj = XjSj 

= (^ +1 -|^ <14) 

This result follows after substituting eq. (13) for X ; and eqs. (4) and (6) 
for ty and tj. Note that Apj is greater than the steady state value igj/A) 
if fj < A and is less than the steady state value if/; > A. 

The expected increase in available pairs next year in the analog entity 
is obtained by summing eq. (14) over all routes, 

n»- = (5f+»(x-2)) 

-!♦'(-£) 

where F is the average fill at next relief for the original entity. The pre- 
vious steady state result [G/A, eq. (12)] is adjusted by a factor that de- 
pends on the change in fill at relief. The adjustment is positive if F < A 
and negative if F > A. 
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4.2.2 Impedance to a change In fill at relief 

The last result [eq. (15)] was derived by assuming an instantaneous 
change in fill at relief from A to fj in all routes in the analog entity at the 
beginning of next year. In practice, we should not expect an instanta- 
neous change. It is not practical or economical to defer a job (by in- 
creasing the fill at next relief) which is already under construction or for 
which material has already been ordered. On the other hand, it is not 
possible to suddenly advance a job (by decreasing the fill at next relief) 
that had been planned for further into the future because major feeder 
relief jobs must be engineered a year or more in advance to allow time 
for job approval, ordering material and construction. 

These considerations can be accounted for in the model by introducing 
the concept of the impedance of the analog entity to a change in fill at 
relief. The effect of this impedance is to cause the fill at relief to be 
changed gradually rather than instantaneously in all routes. In this case, 
the average fill at next relief (F') for the analog entity is somewhere be- 
tween A and F— the average fill at next relief for the original entity. 
Hence F' can be expressed as 

F' - OA + (1 - B)F (16) 

for some 9 between and 1. The parameter 9 is related to the impedance 
of the analog entity to the change in fill at relief. If 9 = 0, the analog entity 
offers zero impedance and the change is instituted instantaneously in 
all routes {F f ■ F). If ■ 1, then the analog entity offers infinite im- 
pedance and the change is never instituted (F' = A). 

The expected available pair increase in the analog entity next year 
based on F' rather than F is 

G ( F'> 

=f +(1 - w ,( 1 -f) (17) 

Note that 9 attenuates the effect of the change in fill at relief on the ex- 
pected increase in available pairs. Equation (17) reduces to our previous 
result [eq. (15)] if the fill at relief is changed instantaneously in all routes 
(9 = 0). An empirical value of 9 will be determined in Section V based 
on the observed behavior of actual entities to a change in fill at relief. 

4.3 Change In average cable size terminated 

In all of the results presented so far, we have assumed that the same 
sized cable is terminated whenever a particular feeder route is relieved. 
In practice (see Ref. 2), changes in many factors can cause future cable 
sizes to be different from those placed in the past. In this case, the new 
cables terminated next year would cause a change in average cable size 
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terminated from its value at the beginning of the year (S). We will now 
derive another adjustment to the model by determining the impact of 
a change in average cable size (AS) on the expected increase in available 
pairs next year. 

The change in average cable size terminated next year can be defined 
in the original entity as 

*S~j}Zxj(cj-Sj) (18) 

where N is the number of routes in the entity, xj is the relief indicator 
function for next year, Sj is the cable size terminated at last relief and 
Cj is the cable size to be terminated at next relief in route;'. We will as- 
sume that a change in average cable size next year is forecast directly for 
the entity. 

If the cable sizes (c, : ; = 1,N) are carried over to the routes of the 
analog entity, then the expected increase in available pairs next year can 
be expressed as 

= ZXj(Sj + Cj -Sj) 

= HXjSj+ Exjicj -Sj) 

= ExjSj + NAS 

The first term is the increase in available pairs in the analog entity if 
there is no change in average cable size next year. The model for its ex- 
pected value is given in eq. (17). If average cable size is forecast to change 
by AS next year, then the increase in available pairs will change by 
NAS. 

The general model is obtained by adding NAS to eq. (17): 

AP(0) - -+ (1 - 6)P ( 1 - -j + AT AS (19) 

This result represents the final adjustment to the original steady state 
model of G/A. The general model can be applied recursively to obtain 
the expected increase in available pairs for any future year. 

4.4 Fill at last relief model 

When defining the analog entity, we assumed that the fill at last relief 
was the same in all routes. A value for this parameter (A) is needed to 
calculate the expected increase in available pairs next year in the analog 
entity [see eq. (19)]. In this section, we derive an estimate of A in terms 
of known aggregate main frame data (see Section 3.2). 

Because A is the fill at last relief in the analog entity, its estimator can 
only depend on parameters which can be observed at the beginning of 
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next year. These include the basic route parameters (wj, Pj, sj : j = 1,N) 
and the corresponding aggregate parameters (N,W,P,S). 

The basic relationship in the analog entity between the route pa- 
rameters (wj, Pj, sj) and A is expressed in eq. (9), 

wj = Apj - Asj(l - r s j s ) (20) 

Both r" and Wj are modeled as random variables and r s j s is assumed to 
be uniformly distributed on the unit interval. 
Summing eq. (20) over all routes in the analog entity yields 
Nw = AP- ANSQ. - R ss ) 

which can be solved for A to obtain 

m NW 

P - NS(l - R ss ) 
In the last two equations, w = (2wj)/N and 

R ss = ^Y.(sj/SK (22) 

are random variables and N, P and S are observed aggregate parameters 

at the beginning of next year. If w is observed to be W/N at the beginning 

of next year, then a random variable (a) can be defined as an estimator 

of A where 

W 

a = (23) 

P-NS(1-R SS ) 

An estimate of A is then the expected value of a. 

The expected value of a does not have a simple form because the 

random variable R ss — which is analogous to the fraction through the 

relief cycle in a single route — is approximately normally distributed. This 

follows from a generalized version of the Central Limit Theorem, called 

Lindebergs Theorem, 5 for uniformly bounded independent random 

variables. In our case, R ss is the sample mean of the independent random 

variables [(sj/S)r] s : j = 1,N] which are uniformly bounded by k = s max /S 

where s max is the largest manufactured cable size. The expected value 

and variance of R ss can be calculated from eq. (22) and the expected 

value (V 2 ) and variance (V12) of rf, 

G(R ss )=jjZ(sj/S)G(r s j s ) 

_ 1 
~2 

vMR ss ) = ^i I>;7S) 2 var(rf ) (24) 
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An estimate of A, A, can be obtained by replacing R ss in eq. (23) by 
its expected value of V2, yielding 

s W 

A = (25) 

P--NS 

2 

This is a reasonable estimator for A if the variance of R ss is small. A 
bound on var(i? ss ) can be derived from the bound on Sj/S. Recall that 
si/S < k, which together with eq. (24) yields 

k 2 
var(R ss ) < 



12N 



Using a value of k = s max /S = 4 which corresponds to a typical average 
cable size terminated S = 900 pairs and the largest manufactured cable 
size s max = 3600 pairs yields var(i? ss ) < 0.003 for a typical area sized 
entity of 400 routes. Hence for large entities, A [eq. (25)] is a reasonable 
estimator of A. 

This estimator will be used for A in the general model [eq. (19)] for 
estimating the available pair increase next year in the analog entity. A 
confidence interval for the model will be derived in the next section for 
the idealized case of an entity of identical routes. The performance of 
the model in actual entities will be analyzed in Section V. 

4.5 Idealized model accuracy 

A confidence interval for the expected increase in available pairs next 
year will be derived for an entity of ^identical routes where the prob- 
ability of relief next year and the cable size terminated are the same in 
all routes. The calculated error bounds for this idealized case are rep- 
resentative of the best possible model performance in actual entities and 
will be compared to observed model errors in Section V. 

If the probability of relief next year and the cable size terminated are 
the same in all routes, then the increase in available pairs next year can 
be represented as a sum of N independent and identically distributed 
binomial random variables, 

The expected value and variance of AP are 

6(AP) = N\s (26) 

var(AP) = Ns 2 \(l - X) (27) 

where X is the probability of relief next year and s is the cable size ter- 
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minated in all routes. If N is large, AP is approximately normally dis- 
tributed by the De Moivre-Laplace Theorem. 5 

A 6-confidence interval for the relative deviation of AP from its ex- 
pected value can be obtained from the probability statement, 

,| AP-£(AP) , x g 
VI £(AP) I / 
which can be expressed in terms of a unit normal random variable, 



<g(AP) \ 



/| AP-<g(AP) 6(AP) \ 

V|(var(AP))^ C (var(AP))W 



If /3 is selected such that 

|AP-£(AP) 



VlCvarCAP)) 1 / 2 V 



then the error bound e can be expressed as 

a (var (AP)) 1 ^ 2 
€ = /? *<AP) 
Substituting eqs. (26) and (27) for £(AP) and var(AP) yields 

l-X\i/2 



«w-'fisr) 



(28) 



This result will be used in the next section as a scale factor to compare 
observed errors in entities with different numbers of routes and average 
probabilities of relief. 

Error bounds e(iV,\) at the 50percent (5= 0.5, = 0.675) and 90 percent 
(5 = 0.9, (8 = 1.645) confidence levels for an entity of identical routes with 
a probability of relief next year of X = 0.2 are given in Table I as a func- 
tion of the number of routes in the entity. These results are represen- 
tative of the best possible model performance in actual entities. 

V. VALIDATION STUDY 

A validation study was performed to determine the statistical accuracy 
of the model when applied in actual entities. This analysis will also 
provide the means of calibrating the model by empirically selecting a 

Table I — Idealized error bounds for available pair increase model 
(X = 0-2) 

Number of Confidence level 
routes 50% 90% 

100 14% 33% 

200 10 23 

400 7 16 

1600 3 8 
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value for 9 based on the observed behavior of actual entities to changes 
in fill at relief. 

5.1 Validation plan 

Route facility charts provide the data required to perform the vali- 
dation study. These charts are maintained for each route in some areas 
to record the past history of assigned and available pairs and to project 
future growth and relief requirements at the main frame. A typical route 
facility chart is shown in Fig. 4. The assigned and available pairs in the 
route at the main frame are counted at year-end and posted on the fa- 
cility chart. An increase in available pairs in the route caused by termi- 
nating a new cable on the main frame is shown in the month that the 
pairs were made available for service. In this case, an estimate of the fill 
at relief can be made by linearly interpolating between the year-end 
counts of assigned pairs made before and after the service date of the 
job. 

The validation study consisted of applying the model in a number of 
sample entities and then analyzing the statistics of the model errors. The 
aggregate main frame data (N,W,P,S,G,F,AS) required by the model 
can be compiled from facility chart data on each of the routes in a sample 
entity. The equations in Sections 3.2 and 4.3 were used with one excep- 
tion: The average fill at next relief (F) was calculated using only routes 
that were relieved next year. Since historical data were used in the study, 
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Fig. 4 — Route facility chart. 
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actual values could be used for G, AS and F rather than forecast values 
which had the desirable effect of eliminating forecast errors from the 
validation study. 

The general model [eq. (19)] was then used to calculate an estimate 
of the increase in available pairs next year in the sample entity for various 
values of between zero and one. An estimated value (A) of the fill at 
last relief (A) required by the model was calculated using eq. (25). The 
actual model error for each value of 6, 

AP(0) - API 



e(6) = 



AP(0) I 
was then scaled using the equation 

where e(N,\) is given by eq. (28) and X is the fraction of routes in the 
sample entity which were actually relieved next year. The scaling is in- 
tended to make errors comparable in sample entities with different 
numbers of routes (N) and probabilities of relief (X). The scale factor 
is arbitrarily defined to be unity when N = 200 and X = 0.2. 

The distribution of scaled model errors over the sample entities was 
determined for each value of 6. The behavior of the 90 percent point and 
the 50 percent point, or median, of these distributions as a function of 
6 was used as the basis for selecting the best value for 6. The 50 percent 
point of the distribution of signed scaled errors was also studied as a 
function of 6 to help detect model biases. Similar studies were performed 
on subsets of the sample entities to determine if the model could be 
calibrated using a universal value of 6. 

5.2 Validation study results 

Facility chart history on available and assigned pairs from over 650 
routes were gathered from three different areas (see Table II). History 
was sought from year-end 1968 to year-end 1975 but this amount of 
history was not available in all routes. 

Sample entities were defined within an area for each year by selecting 
routes whose growth in assigned pairs during the year was at least one 
percent of the total assigned pairs at the beginning of the year. The 
validation study was restricted to only growth routes because the con- 





Table II — Route data 


Area 


Number of routes 


A 
B 
C 


283 
211 
168 
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cepts of timing and sizing relief are most meaningful in a growing 
route. 

This process yielded 21 sample entities (i.e., 7 different years in 3 
different areas). The actual number of routes in each sample entity are 
shown in Table III. Model errors were calculated as described in the 
validation plan (Section 5.1). The 90 and 50 percent points of the dis- 
tribution of scaled errors are plotted versus in Fig. 5. The 50 percent 
point of the distribution of signed scaled errors — which ideally should 
be zero — is also shown to help detect model biases. 

The best value of 6 for estimating the increase in available pairs is 0.6 
which minimizes the 90 percent point of the distribution of scaled errors. 
At this value of 6, the 50 percent point of the distribution is nearly 
minimized, and the 50 percent point of the distribution of signed scaled 
errors is almost zero. Similar analyzes performed for the 7 sample entities 
in each area yielded values of 6 between 0.5 and 0.7. These results indi- 
cate that it is possible to calibrate the model with a universal value of 
= 0.6. 

The anticipated accuracy of the best model (0 = 0.6) when applied in 

Table III — Number of routes in sample entities 











Year 








Area 


"69 


'70 


'71 


'72 


'73 


'74 


'75 


A 


155 


166 


184 


197 


197 


184 


219 


B 


60 


81 


98 


128 


171 


166 


166 


C 


42 


50 


75 


81 


92 


88 


125 



80 
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Fig. 5— Scaled model errors (N = 200, A = 0.2). 
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Table IV — Anticipated error bounds for the available pair 
increase model (0 = 0.6, A = 0.2) 

Number of Confidence level 

routes 50% 90% 

100 19% 48% 
200 13 34 
400 10 24 
1600 5 12 

an entity of N routes with an average probability of relief of X can be 
expressed in terms of the scaled errors as 

e(N,X) 



e(N,\) = e(6 = 0.6) 



e(200,0.2) 



where e(N,\) is given in eq. (28) and e(d = 0.6) is read from Fig. 5 for ei- 
ther the 50 percent or 90 percent confidence level. Anticipated error 
bounds as a function of the number of routes in entities with a probability 
of relief of 0.2 are summarized in Table IV. These results are somewhat 
higher than the corresponding error bounds which were calculated for 
an entity of identical routes (see Table I) because of route-to-route 
variations in the sample entities. 
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